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ABSTRACT 

The first stars forming in minihaloes at redshifts greater than 20 may have been very 
massive, and could have left behind massive black hole (MBH) remnants. In a previous 
paper we investigated the hierarchical merging of these MBHs and their associated 
haloes, using a semi-analytical approach consisting of a hierarchical merger tree al- 
gorithm and explicit prescriptions for the dynamics of merged substructure inside a 
larger host halo following a merger. One of the results was the prediction of a number 
of MBHs orbiting throughout present-day galactic haloes. In addition, we estimated 
the mass-accretion rate of these MBHs, assuming that they retained around them a 
core of material from the original haloes in which they formed. On the basis of these 
estimates, in this paper we determine the bolometric, optical and X-ray luminosity 
functions for accreting MBHs, using thin disk and advection dominated accretion flow 
models. Our predicted MBH X-ray fluxes are then compared with observations of 
ultra- luminous X-ray sources in galaxies. We find that the slope and normalisation 
of the predicted X-ray luminosity functions are similar to those observed, suggesting 
that MBHs could account for some fraction of these sources. 

Key words: galaxies: formation - galaxies: haloes - galaxies: nuclei ~ cosmology: 
theory 



1 INTRODUCTION 

, There is increasing speculation that the first generation of 
stars in the universe may have been extremely massive, and 
that some of these objects could have collapsed directly to 
massive black holes (MBHs) at the end of a brief stellar 
lifetime. If this is indeed the case, than it has interesting 
implications for the formation of the super-massive black 
holes (SMBHs) and unusual X-ray sources observed in the 
local universe. 

Recent semi-analytic (Hutchings et al. 2002; Fuller & 
Couchman 2000; Tegmark et al. 1997) and numerical stud- 
ies (Bromm, Coppi & Larson 2002; Abel, Bryan & Norman 
2000) suggest that the first stars in the universe formed 
inside dense baryonic cores, as they cooled and collapsed 
within dark matter haloes at very high redshift. For a stan- 
dard ACDM cosmology, these minihaloes are estimated to 
have had masses in the range Mmin ^ 10^-10^ Mq , 
and to have collapsed at redshifts Zcoiiapse ~ 20-30 or higher. 
Since these first star-forming clouds contained essentially no 
metals, gas cooling would proceed much more slowly in these 
systems than in present-day molecular clouds. As a result. 
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they may have collapsed smoothly and without fragmenta- 
tion, producing dense cores much more massive than the 
proto-stellar cores observed in star-formation regions today. 
Assuming nuclei within these cores accreted the surrounding 
material efficiently, the result would be a first generation of 
protostars with masses as great as 10"^ Mq (Bromm, Coppi 
& Larson 2002; Omukai & Palla 2001). 

As yet, nothing definite is known about the subsequent 
evolution of such objects. However, their large masses would 
probably result in many of them ending up as MBHs with 
little intervening mass loss - for systems in this mass range, 
gravity is so strong that there is no ejection of material from 
the system in a final supernova bounce (Heger et al. 2002). 
This high-redshift population of MBHs is particularly inter- 
esting, since it might seed the formation of the SMBHs seen 
at the centres of galaxies at the present-day. 

In a previous paper (Islam, Taylor & Silk 2003, paper I 
herafter) we used a semi-analytical model of galaxy forma- 
tion to follow the evolution of MBHs as they merge together 
hierarchically, along with their associated haloes. The code 
we used combines a Monte-Carlo algorithm to generate halo 
merger trees with analytical descriptions for the main dy- 
namical processes - dynamical friction, tidal stripping, and 
tidal heating - that determine the evolution of merged rem- 
nants within a galaxy halo. We introduce seeds into the code 



2 Islam R.R., Taylor J.E. & Silk J. 



Table 1. The mass of seed MBHs, and the height of the peaks in 
the initial density field in which they formed, for the four models 
considered. The collapse redshift ^collapse 's the epoch when peaJjs 
of height i^pfc first cross the cooling threshold (see paper I) . 



Model 




peaJj height Vp}^ 


^collapse 


A 


260 Mq 


3.0 


24.6 


B 


1300 M0 


3.0 


24.6 


C 


260 Mq 


2.5 


19.8 


D 


260 M© 


3.5 


29.4 



by assuming that in each minihalo forming before a redshift 
•^collapse, a single MBH forms as the end result of primordial 
star formation. 

For our computations, we considered four different sets 
of values for the parameters v^k and M,^seed, which fix the 
abundance and the mass of the seeds respectively. These 
values are summarised in table 1. Our choice of a seed MBH 
mass of 260 Mq is motivated by the result of Heger (2002) 
that massive stars above this mass will not experience a 
supernova at the end of their lives, but instead will collapse 
directly to a MBH of essentially the same mass. 

In paper I we investigated the abundance of MBHs in 
present-day galaxies as a result of this process of hierarchical 
merging and dynamical evolution. This is shown in figure 1. 
We now look at how this model could be tested observation- 
ally in nearby galaxies. As the MBHs orbit they will emit 
in the optical, UV and X-ray due to accretion of matter in 
their direct vicinity. This matter can be either the interstel- 
lar medium (ISM) of the host, or remnant material of the 
original satellite that remains bound to the MBH. While op- 
tical and UV emissions could in principle arise from a variety 
of sources and be subject to a range of absorption and dis- 
tortion mechanisms, X-ray emissions are more uniquely as- 
sociated with accreting black holes. In addition, hard X-rays 
are much less subject to absorption. Observations of lumi- 
nous and ultra-luminous X-ray point sources associated with 
galactic haloes in the local universe might therefore provide 
constraints on the existence and abundance of a population 
of MBHs. Optical and UV observations could be used as a 
follow up, to confirm the CEise for detection of a MBH. 

In section 2 we outline the accretion scenarios and mech- 
anisms we considered. These include thin disk and advec- 
tion dominated accretion flows in particular. In section 3 
we present the resulting optical and X-ray luminosity func- 
tions for MBHs orbiting in present-day haloes. This is then 
compared with observations of ultraluminous X-ray sources 
(ULXs) in section 4. We sum up our results and conclude in 
section 5. 



2 ACCRETION PROCESSES 
2.1 Accretion rates 

We assume that accretion onto MBHs orbiting in present- 
day haloes is given by the Bondi-Hoyle mass-accretion rate 
(Bondi & Hoyle 1944; Bondi 1952) 

dM 



dt 
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Figure 1. Abundance of all MBHs in the halo for models C, A 
and D (top to bottom panels) for different final halo masses aver- 
aged over 30 trees with error bars corresponding to the standard 
deviation. 



Here Ca is the sound speed in the gas and pg its density - 
both far from the MBH. v, is the velocity of the MBH and 

race the accretion radius 
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where we have used /3s = v,/cs (Chisholm, Dodelson & Kolb 
2003). If the MBH accretes adiabatically from a gas of pure 
hydrogen this is 



M = 8.77 X 10" 
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If a constant fraction rj of the accreted mass is radiated away, 
the bolometric accretion luminosity of the MBH is then 



Lboi = 



dM 2 



Pa 



/ M. 

\ IOOM0 J \ iU-^" g cm 
T)"'(l + /3?r/%rgs-^ 



(5) 



V 10 km s" 

where c is the speed of light. 

In what follows we assume that while the nature of the 
accretion process is somewhat uncertain, the overall mass- 
accretion rate is essentially determined by the Boridi-Hoylc 
formula. We neglect the possibility that the mass-accretion 
rate is modified e.g. by a non-negligible mass of an accretion 
disk that may form around the accreting MBHs. 

As the MBHs orbit through the host halo they accrete 
matter from the host ISM. In paper I we have shown that 
ISM accretion rates are too small to lead to significant accre- 
tion signatures even in regions with relatively large amounts 
of gas, such as the disk or the galactic centre. Much higher 
rates are possible if MBHs accrete from a residual core of 
material from the satellite they were originally associated 
with. Such a core, consisting of gas and possibly stars, could 
remain bound even if the original satellite has lost all but 
a few percent of its initial mass. While insignificant for the 
overall dynamics of the satellite, the amount of baryonic 
matter contained in this core could still contribute signifi- 
cantly to the accretion onto a MBH that is present at its 
centre and so boost its accretion luminosity, potentially by 
orders of magnitude. In this case MBH accretion is indepen- 
dent of local ISM density and of the velocity of the MBH 
relative to the host halo, as the 'fuel supply' for accretion 
travels with the black hole. 

To get a sense of the possible mass-accretion rates al- 
lowed by the mechanism, we need to estimate the gas den- 
sity in the residual baryonic core. We will assume that all 
the baryons in the original system cool and collapse down to 
a radius of order ~ Q.lr^ir (where r^ir is the virial radius 
of the original system), at which point thermal pressure or 
angular momentum halts the collapse (equivalently we could 
assume that only a fraction of the baxyons collapse, but that 
rb is slightly smaller). If all baryons are in the form of gas 
and the baryon fraction in the satellite is cosmological, i.e. 
Mb = iyib/^m) Mvir, then the mean gas density is 

3Mb 1 Qb SM^ir „„„ 102^6 ^^ir 

Pa = 7—3 = 777^ 7^ « 2.39x10 — -3— (6) 



A-nrl 0.001 D.m ^-Krl^^ 



Substituting this into equation (5) and using Cs « 10 km s~^ 
for ISM at 10* Kelvin we can determine the mass-accretion 
rate 



6.25 X 10" 



(Mmbh'^ 




f M^ir \ 


\ Mq J 




[w^MqJ 



\ —3 

kpc / VlOkms-i/ 



(7) 



and thus the luminosity. Here m is the mass-accretion rate 
in units of the Eddington accretion rate 
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2.2 Accretion mechanisms 

As material gets accreted onto a Black hole, a fraction of 
its gravitational potential energy released as radiation. The 
amount and spectrum of energy released in such a way de- 
pends on the accretion model, that is a prescription for the 
physical processes that facilitate the dissipation of energy 
as material moves towards the accreting MBH. We consider 
two of these. 



2.2.1 Thin disk accretion 

Shakura & Sunyaev (1973) proposed a model in which black 
holes (BHs) accrete from a geometrically thin but optically 

thick accretion disk, which is considered to achieve the high- 
est radiative efficiencies, typically with 'qtd ^ 0.1. Because 
the disk is optically thick, radiation is emitted locally with 
a thermal spectrum and black body temperature 

nr) = {^'i=r (9) 



/ 3GM.M 

where M and r are the mass-accretion rate and radial dis- 
tance from the MBH respectively. We can convert these into 
dimensionless quantities and express the accretion rate in 
units of the Eddington accretion rate equation (cf. equation 
8) and the radial distance from the MBH in units of the 
MBH Schwarzschild radius, r. 
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The corresponding spectrum integrated over the whole 
disk is given by (Frank, King, & Raine 1985) 



fdf 



/kT(f) _ I 



(11) 



where ri^o = 3 r. is the last stable orbit radius of the MBH 
and ^0 = 4 TT is the area of a spherical shell at unit distance 
from the MBH. However, in what follows we treat ^0 more 
generally as a normalisation constant. All our computations 
were carried out with rtd = lOOOr, and ^0 was consequently 
renormalised by setting ^ ( dLu/ Ai>)Av = Lboi. The result- 
ing spectra are shown in the panels on the right in figure 
2. 

In general the accretion flow near the BH may become 
optically thin and absorption of photons becomes less im- 
portant than scattering, which results in an increase of the 
mean energy of the photons radiated away. The emergent 
flux is thus non-thermal and has a higher temperature. Con- 
sequently in the inner parts of the disk near the MBH, 
equations 9 and 11 are not strictly valid anymore. Another 
prescription to account for the higher frequency photons is 
needed. 

This high-frequency emission also has important obser- 
vational consequences. MBHs are only clearly identified as 
BHs via an X-ray signature. However, from equation 9 the 
peak frequency of the photons from a standard thin disk 
scales T{ruo) oc A^y**. For the MBHs any resul- 

tant fiux in the X-ray range is only significant at energies 
below about 1 keV. A standard thin disk alone would there- 
fore not provide the X-ray fiux by which to identify MBHs. 

One could tackle this problem phenomenologically by 
adding a power-law spectrum extending into the X-ray 
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range. A power law X-ray spectrum is well established for 
active galactic nuclei (AGN), which are driven by SMBHs, 
and is also found in recent observations of ultra-luminous X- 
ray sources in nearby galaxies, which could well be accreting 
MBHs (Colbert & Mushotzky 1999; Wang 2002). However, 
instead we consider a different model that has recently been 
developed particularly with a view to explaining the emis- 
sion spectra from accreting SMBH in low-luminosity AGN. 
This is an advection dominated accretion flow. 



2.2.2 Advection dominated accretion flows (ADAFs) 

In the ADAF model, once the accretion rate drops below 
some critical rate which we define below, accretion 

switches to an 'advection dominated accretion flow' (ADAF) 

(Ichimaru 1977; Rcos ot al. 1982; Narayan & Yi 1994; Man- 
moto, Mineshige & Kusunose 1997). Because of the low ac- 
cretion rates the density of the accreted gas is low, too. The 
gas, not being able to cool, efficiently stores thermal energy 
instead of radiating it away and carries it into the BH. On 
the other hand, because the accretion flow is optically thin 
the spectrum of the emergent radiation is modifled by scat- 
tering processes. The main processes are synchrotron emis- 
sion, brcmsstrahlung and comptonisation. For more details 
the reader is referred to the review by Narayan, Mahadevan 
& Quataert (1998) and references therein. 

The spectrum emitted from a BH accreting in this way 
depends on a range of parameters such as the mass-accretion 
rate , the viscosity pararnctcr and the relative contribution 
of magnetic fields. Typically the spectrum has to be deter- 
mined numerically, but can be approximated by analytic 
scaling relations (Mahadevan 1997). Using these, the total 
radiative efficiency for an ADAF, rjadaf can be expressed in 
terms of that of a thin disk, rjtd as 



Vadaf = 0.2rit, 



9{0e) 



(12) 



for a typical set of parameters as given in Mahadevan (1997). 
Here m is the mass- accretion rate in terms of the Edding- 
ton accretion rate, a is the viscosity parameter (a ~ 0.1 to 
0.5), 9e = fcs Ts/{me c?) is the electron kinetic energy in 
units of their rest mass energy. g(6e) ~ 11.5(Te/10®K)~"'^'^ 
is a measure of whether electrons in the flow move at rel- 
ativistic speeds or not. Typically 5 ~ 1 to 10 for electron 
temperatures 2 x 10^ <Te < 10^° K. We can determine an 
equilibrium electron temperature which is approximately 
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Generally the effect of ADAFs is to reduce the overall 
radiative efficiency compared to that from thin disk accre- 
tion, although individual regions of the accompanying spec- 
trum, especially the X-ray band, tend to be much stronger 
than for a thin disk. 

As mentioned above the thin disk and ADAF regime 
can be dehneated by the critical mass accretion rate rhcrit 
which is given by 



merit ^ 0.28a 



(14) 



For rates above merit the ADAF mechanism breaks down 
and in this regime we will work with the thin disk model. 



For low values of a 0.1 accretion may switch to a 'con- 
vection dominated accretion flow' (CDAF) (Ball, Narayan 
& Quataert 2001), which is characterised by much stronger 
suppression of mass accretion. We will neglect this possi- 
bility here as we work with a = 0.3 in the computation of 
spectra and luminosity functions for our model. This is also 
the fiducial value used in the work of Mahadevan (1997). We 
will briefly consider changes in a below. 

In paper I we have seen that MBHs do not accrete at 
more than 1/10 of the Eddington rate with most accreting 
at less than 1/100. Thus the ADAF model can be applied 
to most MBHs we are dealing with. 

The three dominant contributions to the continuum 
spectrum of an ADAF are those from synchrotron and 
brcmsstrahlung emission as well as comptonisation, which 
can all be separately approximated by scaling laws (Mahade- 
van 1997). In the radio-submillimeter range synchrotron 
emission dominates and the luminosity (per unit frequency) 
scales as 



M, 



a \ -4/s / T, 



2/5 



M, 



up to a peak frequency Vp given by 



109 K, 

s"' Hz"' (15) 

"1/2 
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X [-^ ) m ' Uz (16) 



Photons from synchrotron emission can get comptonised, 
which raises their energy up to a limit given by the elec- 
tron temperature fcs Tg. Since most synchrotron photons 
are emitted at v = fp we can make the simplifying assump- 
tion that this is the initial frequency for all photons to be 
comptonised. The comptonisation spectrum is then simply 



dL. 



compt 



dv 



(17) 



where the power law index is primarily a function of electron 
temperature and mass-accretion rate 
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(18) 



There is a constant contribution from bremsstrahlung 
emission that also tails off exponentially at the max- 
imum frequency implied by the electron temperature. 
Brcmsstrahlung only becomes prominent in the UV/soft X- 
ray range if Oc > 1, such that power from Compton emission 
effectively decreases with increasing frequency 
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If < 1, F{e) takes the form 
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whereas for the case 6, > I, 

Fie,) = (^^^ [ ln(1.123ee + 0.48) + 1.5] 
+2.3ee{ lnl.l236'e + 1.28) 
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The resulting spectrum is shown in the left panels of figure 
2 for a range of MBH masses and accretion rates ranging 
from m — rhcrit to m = 0.03 rhcrit, where Thcrit = 0.0252 
for a = 0.3 which we use for all our computations. 

We see that only ADAF accreting MBHs display sig- 
nificant X-ray emission. It is therefore mainly MBHs that 
could be identified as such by their X-ray signature. Con- 
versely, while optical and UV observations may be used as 
follow-ups, the sources that are most luminous in this range 
are accreting from a thin disk and so will not have any as- 
sociated strong X-ray emission. 



X-ray hardness ratios 

For our predictions we have computed the X-ray emissions 
for the energy range between 0.5 < h v < 2.0 keV and 
2.0 < h u < 10.0 keV , which we call the soft and hard X- 
ray bands respectively. Both bands lie in the Bremsstrahlung 
part of the ADAF spectrum where dL/ di/ = 1, i.e. the pho- 
ton counts from both bands are the same; their hardness 
ratio is equal to 1. Unfortunately, observations from X-ray 
binaries (XRBs) can yield similar values, so it will not be 
easy to distinguish between MBHs and XRBs on the basis 
of this alone. The combined thin disk -I- ADAF model that 
we used, means that we can also have sources with large 
soft X-ray luminosities and very small hardness ratios. This 
is the case particularly for seed MBHs for which thin disk 
accretion reaches furthest into the soft X-ray regime, and 
there actually dominates over the respective ADAF contri- 
bution (cf. figure 2). In section 3.2 we see that this leads 
to a distinctive 'bump' in the soft band luminosity func- 
tion. These sources should be identifiable by their very low 
hardness ratio. 



3 PREDICTED LUMINOSITY FUNCTIONS 

Using the spectral models outlined above we have computed 
optical and X-ray luminosities of accreting MBHs in present- 
day haloes. For comparison we show in figure 3 the bolomet- 
ric luminosity functions, obtained in paper I. 



3.1 Optical emission 

For the MBH masses we are considering here, the optical 
region lies within the Rayleigh- Jeans part of the extended 
black body spectrum that characterises the corresponding 
thin disk emission. In this part, kBTljtd) hf <^ kBT{riso) 
where the spectrum can be approximated by a power law 




42 38 40 42 
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Figure 3. Bolometric luminosity of MBHs accreting from bary- 
onic cores for models C, A and D (left to right panels). Shown 
are the sources whose line of sight falls within some projected 
distance from the host centre. 
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Figure 4. Average B-band luminosity functions for MBHs in 
haloes of various masses, for two different values of a (two sets 
of lines). The dashed lines show the contribution from ADAF 
accretion. Thin disk emission dominates at large luminosities. 
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Figure 2. Emission spectra from ADAF (left) and thin disk (riglit) accreting MBHs in the mass range between 2.6 X 10^ to 2.6 X 10^ 
M0 . In each panel the curves from top to bottom correspond to mass-accretion rates rn = {1, 0.3, 0.1, 0.03} merit, where merit = 0.0252 
for a = 0.3. 



and the constant A\ depends on the MBH mass and accre- 
tion rate as 



^1 = 2.57 Ao h-"'"' c-^ T{ruof"' 

^ 1.64 X 10" Ao m-l'^ M^^ erg s~^ Hz" 



(23) 



In the ADAF regime it is primarily synchrotron emis- 
sion and comptonisation that contribute to the emission in 
the optical/UV range, with the corresponding scalings given 
in equations 15 and 17. 

Figure 4 shows the B-band luminosity function. The 
dashed line represents the contribution from ADAF emission 
with the rest, primarily at large luminosities, coming from 
thin disk emission. Apart from our fiducial value for the 
viscosity parameter, a = 0.3 (thick lines), we have also given 
results for a = 0.5 (thin Unes). While the results for ADAF 
emission are typically not too sensitive to changes in ADAF 
parameters, results are affected more strongly by changes in 
a, but only through its direct impact on the value of merit 
(cf. equation 14). As a result the most important effect of 



increasing a is to extend the number of ADAF emitting 
sources to larger luminosities. Since a thin disk radiates more 
efficiently than an ADAF this results in a lower abundance 
of sources at very large magnitudes in the optical range. For 
smaller values of a we get a correspondingly larger number 
of thin disk accreting sources. However, we found that for 
a value of a — 0.1, for instance, the resulting luminosity 
function in the B-band displays a significant 'dip' which is 
also present to some degree for a = 0.3 in the case of the 
most massive 1.6 x 10^^ Mq halo in figure 4. We will deal 
with the results for the X-ray range in the next section, but 
note at this point that a higher (lower) value of a would 
result in a larger (lower) abundance of sources in the X-ray 
regime, emission in which is mainly produced in ADAFs. 

In figure 5 we show the B and V band luminosity func- 
tion of MBHs within some projected distance from the host 
centre normal to the line of sight. Here we have included 
all MBHs within the respective host virial radius along the 
lines of sight passing within some projected distance from 
the centre. We will subsequently refer to this as the projected 
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luminosity function. The projected luminosity function can 
more easily be compared with observations, which typically 
refer to sources within galaxies. In each case the function is 
plotted for MBHs sources within Rproj =0.1, 0.03, 0.01 Rvir 
from the host centre, whore the latter corresponds roughly to 
the outer radial extent of the gala^xies in the haloes. Across 
all halo masses we find that a fraction of about 30, 9 and 0.9 
% respectively of all haJo MBHs are located within the pro- 
jected distances listed above. Since there is no discernible 
difference in the relative distribution of MBH masses at dif- 
ferent radii (paper I), this means that for any one MBH in 
a given mass, accretion rate or luminosity range and within 
Rproj ~ O.lRvir, there are about two more in the halo. We 
have marked with bold black arrows the luminosity range 
below which ADAF emission dominates. 

It is important to note that these luminosity functions 
only provide an upper limit on what could actually be de- 
tected. We have here neglected the effect of absorption which 
is expected to significantly reduce the flux of optical and UV 
photons due to dust and neutral hydrogen along the line of 
sight. The latter affects particularly the light from sources 
whose line of sight passes closely to the galactic centre, such 
as those with Rproj < O.OlRvir and located on the far side 
of the halo. We would expect this effect to be less important 
in gas-poor ellipticals. 

Even if absorption is accounted for, the large luminos- 
ity of the accreting MBHs by itself may not be sufficient to 
identify them as such. As they appear as point sources, they 
could also be XRBs or background AGN. For the largest 
accretion rates, m ~ 0.1 (paper I), at which MBHs accrete 
from a thin disk in our model, it is possible that emission oc- 
curs in the form of (mini) jets which may be observable even 
though the optical and UV output from across the accretion 
disk may not. 



3.2 X-rays 

In our model X-rays axe primarily emitted by ADAF ac- 
creting MBHs. Emission in the X-ray range is dominated by 
comptonisation or bremsstrahlung depending on the Comp- 
ton spectral power law index, ac (eq. 18). The two contri- 
butions are determined by the scaling relations equations 17 
and 19. In figure 6 we present the projected X-ray luminos- 
ity function, both in the 0.5 - 2.0 keV and 2.0 - 10.0 keV 
energy range. Again we have ignored the issue of absorption 
here, which can be important primarily for the softer 0.5 - 
2.0 keV X-ray band. 

X-ray emissions are the primary signature to look for 
in observations, while optical and UV detections are prob- 
ably useful only for follow-up, due to the range of poten- 
tial sources and the problem of absorption. Another sec- 
ondary diagnostic would be to see whether any correlations 
exist between different parts of the emitted spectrum. How- 
ever, this would be extremely dependent on the accretion 
model. For instance, in the ADAF model that we used, the 
shape of the spectrum near the exponential cut-off in the 
X-ray regime and in the radio/submillimeter range is essen- 
tially unaffected by the mass and accretion rate of the MBH 
(cf. figure 2). If a MBH accretes via an ADAF we would 
therefore expect a correlation of the radio and X-ray lumi- 
nosities. 



4 COMPARISON WITH OBSERVATIONS 

We now compare our predictions with observations. As we 
have already noted. X-ray emissions will be the most impor- 
tant trace to look for. We will therefore confine our compar- 
ison to X-ray observational data only and further more to 
sources which have been identified as belonging to a galaxy. 

In our predictions we have generally not taken into ac- 
count absorption of MBH emission either in the originating 
galaxy as well as ours. To minimise the former one would 
want to look for sources in external galaxies that we see 
face-on. This is especially important for gas-rich spirals. To 
account for the extinction from the IGM and the Milky- 
Way's ISM all-sky extinction maps can be used (Schlegel, 
Finkbeiner & Davis 1998). These allow a determination of 
the average hydrogen column density between us and the 
target galaxy. 

4.1 Observations of ultra-luminous X-ray sources 
(ULXs) 

For the baryonic core densities and accretion mechanisms 
considered above, our model predicts the presence of ultra- 
luminous X-ray sources (ULXs) in galaxies and their haloes. 
By ULX we refer to any compact X-ray source with Lx ^ 
10^* erg s~^. This cutoff is interesting because it corresponds 
roughly to the Eddington limit for isotropic accretion onto 
stellar-mass black holes. Early observations revealed that 
some galaxies contain these very bright sources outside 
their nucleus (Fabbiano 1989). Subsequently, a number of 
other studies discovered ULXs in data from the ROSAT 
and ASCA satellite missions (Colbert & Mushotzky 1999; 
Makishima et al. 2000; Roberts & Warwick 2000; Colbert 
& Ptak 2002). With the advent of the Chandra and XMM 
satellite missions, the number of ULX detections has soared. 
The superior resolution and sensitivity of these satellites 
has greatly increased the distance at which ULXs can be 
detected, and thus the number of galaxies in which their 
properties can be studied. 

While 'normal' stellar remnants may easily account for 
the fainter ULXs, sources which emit more than a few 
times lO'^^'erg require black holes more massive than 
those produced by conventional stellar populations, assum- 
ing the emission is Eddington limited and occurs isotropi- 
cally. If ULXs do not emit isotropically, however, but pro- 
duce beamed emission that happens to be aligned with 
the line-of-sight, then stellar mass BHs can be the source 
(King et al. 2001; Zezas & Fabbiano 2002). A departure 
from isotropy implies that the Eddington luminosity is no 
longer a limiting factor, either, as radiation pressure need 
not act along the same axis as the accretion flow. High mass 
X-ray binaries (HMXRBs) have been suggested as candi- 
dates for anisotropic ULXs. Within galaxies, ULXs appear 
to be found predominantly in star-forming regions (Kilgard 
et al. 2002; Zezas & Fabbiano 2002), i.e. they seem to be 
associated with starburst galaxies or the spiral arms of disk 
galaxies. For a number of these galaxies it was found that 
the ratios of ULX number to massive star formation rates 
are similar (Smith & Wilson 2003). Thus many ULXs may 
be young HMXRBs, which would naturally be found in or 
close to star-forming regions. 

There is some evidence, however, to suggest that not 
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Figure 5. Optical luminosity function for MBHs accreting from baryonic cores within various projected distances from the host centre. 
Results are shown for models C, A and D (three sets of panels, left to right) in the B and V-band (left-hand and right-hand figures). 
The bold arrows delineate where luminosities are dominated by thin disk emission and ADAF emission respectively. 
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Figure 6. X-ray luminosity function of MBHs with (linc-of-sight within) various projected distances from the host centre (see 
Results are shown for models C, A and D (three sets of panels, left to right) and X-rays in the energy range from 0.5 - 2.0 keV 
figure) and 2.0 - 10.0 keV (right-hand figure). 



also text), 
(left-hand 



all ULXs fall into this category. Multi-colour disk (MCD) 
fits to ULX spectra indicate that the maximum colour tem- 
perature associated with some sources is significantly lower 
than would be expected for a stellar mass BH system. If 
the maximum temperature arises from emission close to the 
inner radius of the accretion disk, and we assume this to 
be the BH's last stable orbit, a lower temperature implies 
a larger stable orbit and thus a more massive BH (Wang 
2002). Furthermore, ULXs have been detected in some gas- 
poor ellipticals (Jeltema & Canizares 2003), and Colbert & 
Ptak (2002) find that the number of ULXs per galaxy is ac- 
tually higher for ellipticals than for non-ellipticals. Thus the 
connection between ULXs and star formation established 
for spiral galaxies may not be universally valid. It seems in- 
teresting, therefore, to consider alternative explanations for 
these sources, such as isotropically emitting MBHs. 

If seed MBHs do account for some of the observed 
ULXs, there are a few remaining puzzles to work out. It 
is not clear whether our model should predict more ULXs 



around spirals or ellipticals, and the morphological informa- 
tion in the model is probably too crude to test this defini- 
tively. Since ellipticals are generally older systems, however, 
it is not implausible that they would have a large num- 
ber of ancient satellites in orbit around them. A number 
of ULXs seem to be located in globular clusters (Jeltema & 
Canizares 2003), suggesting binaries are the corresponding 
ULX sources. In the context of our model, however, some of 
these globular clusters might be precisely the baryonic cores 
we have considered previously. 

Ultimately, detailed spectral information should allow 
us to distinguish between binaries and MBHs as ULX candi- 
dates. Reasonable spectral information is available for some 
systems, although fits to such data still assume model spec- 
tra which should be verified independently. For the more 
distant sources, the detected count rates are so low that 
spectral fits cannot be carried out. Instead a hardness ratio, 
the ratio of the number counts in a soft and a neighbouring 
hard X-ray band, can be used as a crude diagnostic. Uncer- 
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tainties in the correction for absorption in soft X-ray bands 
complicate the interpretation of these results, however. Here 
we will avoid these complications, using spectral information 
only to the extent that it allows X-ray band luminosities to 
be derived from raw photon counts. Wo will not use spectral 
fits to make any inference about the ULX source object. 

In summary, while binaries may explain the majority 
of ULXs, especially in star-forming galaxies, MBHs are not 
ruled out as a possible source of some of these objects. In 
what follows, wc will show how the frcc-floating MBHs pre- 
dicted by our model could account some fraction of all ULXs. 
The exact fraction is highly model dependent; for the param- 
eters assumed above, we will show it can be considerable. 

4.2 Comparing Luminosity Functions 

Table 2 gives a summary of the observations of point sources 
in individual galaxies that we are considering. Their inferred 
mcisses roughly place them in a category with the 1.6 x 10^^ 
Mq halo in our model. Observed sources are considered 
roughly within the light radius of the galaxy. Where ob- 
servations do cover regions significantly outside this radius 
it is typically an area to one side of a galaxy that happens 
to be still covered by the detector. Interestingly, a number 
of ULXs are detected in this area as well. Estimates of the 
expected number of background objects, such as AGN, su- 
pernovae (SN), etc. show that they can account for most of 
the ULXs observed outside the light radius, although the er- 
ror margins are large. A number of sources outside the bulge 
have projected locations that appear to be in the disk. In 
this case background sources cannot account for the rela- 
tively large number of sources that seem to lie in the disk. 
It is assumed on this basis that these sources are actually 
located in the disk, where their number can be explained 
more naturally. For sources with apparent locations in the 
disk as well as outside the light radius our model offers the 
alternative explanation that some of these sources could also 
be located throughout the halo but inside a column along 
the line of sight that projects onto the disk. 

Given the uncertainties and free parameters in our 
model, wo cannot place a strong constraint on MBHs from 
ULX observations alone. Nonetheless, it is interesting to 
make a quantitative comparison between the observations 
and our model predictions, particularly since the predictions 
scale in a straightforward way as the number or mass of the 
seeds, or the density of the baryonic cores, varies. Since our 
models often predict one or fewer luminous sources per sys- 
tem, the comparison with observations will be statistical; 
this is all the more true since we predict a large halo-to- 
halo scatter in the number of MBHs, and thus in the X-ray 
luminosity functions. In general, even if MBHs accounted 
for all ULXs, we would expect a variation in the luminosity 
functions of individual galaxies that was comparable to the 
halo-to-halo scatter in our results. 

We first look at how our results compare with obser- 
vations of sources in individual galaxies. Figure 7 shows the 
cumulative X-ray luminosity function from observations and 
those that we predict for a 1.6 x 10^^ Mq halo, which would 
be expected to host galaxies of the mass wc arc compar- 
ing with. For our predictions we have chosen two different 
X-ray bands as shown in the figure to match those from 
the observations more closely. Our predictions are consistent 



with the data from the individual galaxies if we only con- 
sider the ADAF emissions. Inclusion of the thin disk compo- 
nent leads to a much flatter slope and an over prediction of 
sources at the high luminosity end. The flat slope essentially 
bridges the gap in the differential luminosity function be- 
tween ADAF accreting MBHs and a small number of MBHs 
accreting from a thin disk as shown in figure 6. Whether the 
thin disk is included or not, MBHs cannot account for the 
large number of sources at low luminosities, which is not 
a problem as we expect XRBs to certainly dominate this 
regime. 

For luminosities larger than Lx ~ 10^* erg s~^, the pre- 
dicted average numbers of objects are near or below 1 and 
wc have to compare with observations of a statistical sam- 
ple of galaxies. The two data sots wc use for this arc that of 
Roberts ^Warwick (2000, hereafter RWOO) and Colbert & 
Ptak (2002, hereafter CP02) . Both sets of observations have 
been obtained using the ROSAT High Resolution Imager 
(HRI) instrument. In comparing with these observations, we 
quote the error on the mean (and not the standard deviation 
of the halo-to-halo scatter) for our predictions. The two sets 
of observations vary fundamentally in that CP02 only con- 
sider galaxies that do contain ULXs. They find a total of 87 
sources in 54 galax;ies of which 15 are ellipticals. They note 
that in their sample the number of ULXs per galaxy is larger 
for ellipticals, which account for a total of 37 ULXs. This 
corresponds to 2.3 ULXs per galaxy for ellipticals compared 
to 1.3 per galaxy in the case of non-ellipticals. As the CP02 
sample only counts galaxies that contain ULXs the derived 
luminosity and spatial distributions will have a systematic 
offset towards higher abundances. 

The differential and cumulative luminosity functions are 
shown in figure 8, for the observations and for our expected 
source population for the galaxy-sized 1.6 x lO^^M© halo in 
model A. We note that the luminosity bands used in the 
two data sets are different; while RWOO quote results for 
the ROSAT 0.1 - 2.4 kcV band, CP02 give their data for 
the 2.0 - 10.0 keV range. To match these we have slightly 
changed our soft bands from the one used above (0.5 - 2.0 
keV). Furthermore, we have only taken into account sources 
at radial distances less than 10 % of the halo virial radius, 
which corresponds to about i? < 30 kpc for this halo. This 
is also of the order of the limiting radii below which sources 
in the CP02 sample are considered as belonging to a galaxy. 

Between 3x 10^* and 10*° erg s^^, our predictions in the 
hard band and the soft band without the thin disk contribu- 
tion match the slope and normalisation of RWOO. A power- 
law best fit yields a logarithmic slope of — 0.83±0.03,— 0.78± 
0.01 and -1.08±0.01 for our prediction, RWOO and CP02 re- 
spectively. For lower luminosities the down-turn (tailing off) 
in the differential and cumulative distribution of the RWOO 
data is probably due to incompleteness. Even if no MBHs 
were present there would still be enough XRBs to match 
or exceed the number of objects at lower energies. Including 
the thin disk in the soft band leads to the black curve on the 
right, which is ruled out even by the CP02 data which, as 
we mentioned above, is biased towards higher abundances. 

Aside from the uncertainty in our model parameters, 
there arc also several physical processes that may reduce 
the luminosity of the brightest sources. The relatively large 
accretion rates in MBHs accreting from thin disks could ex- 
haust the baryonic core, if the latter was not large enough 
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Table 2. Summary of X-ray observations. The first four galaxies are spirals; NGC 720 is an elliptical. All observations were carried out 
with the ACIS instrument on Chandra, except for the observation of M31 which was observed with Chandra's HRC instrument. 





Energy range 


Region surveyed 


# ULXs 


Reference 


M31 (S) 


0.1-10 keV 


Rproj < 1.13kpc 


3 


(Kaaret 2002) 


M81 (S) 


0.3 - 8 keV 


Rproj > D25 


2 


(Swartz et al. 2003) 






disk 


10 








bulge 


5 






0.2 - 8 keV 


disk {Rproj < 8.5kpc) 


5 


(Tennant et al. 2001) 






bulge 


3 




NGC 6964 (S) 


0.5 - 5 keV 


Rproj < 13kpc 


15 


(Holt et al. 2003) 


NGC 1068 (S) 


0.4 - 5 keV 


Rproj ^ 


40 


(Smith & Wilson 2003) 


NGC 720 (E) 


0.3 - 10 keV 


Rproj < 25kpc 


41 


(Jeltema & Canizares 2003) 




[erg s"'] 



Figure 7. Cumulative X-ray luminosity functions for five galaxies observed by Chandra (see text), compared with our model predictions. 
The solid black line with error bars shows our predicted cumulative luminosity functions for model C in the 0.1-10 keV band, while the 
dashed line immediately below it shows the same for model A. The two lines below this are for models C and A in the 0.3-10 keV band. 
In each case only those sources with projected radii < 30 kpc from the galactic centre are included. 



to start with. Past merger activity of MBH hosts might also 
have triggered periods of mini-quasar activity that has lead 
to a faster depletion of gas in the baryonic core. Alternately, 
if the accretion rates were slightly lower than assumed, by a 
factor 5 or so, then most systems would fall below the crit- 
ical accretion threshold (equation 14), and would accrete 
in the ADAF mode. Abandoning the thin disk contribution 



would also remove the brightest sources in the B and V band 
luminosity functions. 

In figure 9 we show the number of sources vs projected 
distance from the centre. The predictions we show are for 
model A and the same parameters as in figure 8. We did 
omit the results for the soft band without the thin disk, as 
this curve is essentially the same as the one for the hard 
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Figure 8. DifTercntial (top panel) and cumulative (bottom) luminosity functions for predictions from model A (1.6 X IO^^Mq halo) 
and observations. Note in particular the similar slope for luminosities larger than lO^^erg . Only sources at radial distances larger 
than R > 0.01 R^ir ~ 3 kpc were considered. The predicted values only account for sources at projected distances less than Rproj < 
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Figure 9. Projected radial distribution of sources with Lx > lO^^erg s~^, as predicted from model A (1.6 X IO^^Mq halo), compared 
with observations. Only sources at radial distances larger than R > 0.01 R^ir ~ 3 kpc were considered. The predicted values only account 
for sources at projected distances less than Rproj < 0.1 Rvi,- ~ 30 kpc. 



band. The slopes of the spatial distribution of the latter - 
and by extension also that of the soft band without thin 
disk - also match that of the observations reasonably well. 
However, agreement gets increasingly worse for the RWOO 
sample towards small host distances where its slope gets 
considerably steeper. This is probably because RWOO have 
imposed a limiting minimum angular distance below which 
they did not consider any sources to avoid confusion with a 
central source. Just as we found no trend in the distribution 
of MBH masses with radius in paper I, we do not find any 
indication of 'luminosity segregation' as a function of radius 
here. 

Apart from the systematic bias in the CP02 data, the 
normalisation is also affected by varying degrees of com- 
pleteness in the two samples. For instance the CP02 set 
features no source with Lx < 10^" erg s"\ and the RWOO 
data are restricted to much smaller radii. Both CP02 and 
RWOO maintain that contamination with background AGN 
is rather small and accounts for no more than about 11-13 
per cent of the sample (RWOO). In the near future, larger 
compilations of data should be available from XMM and 
Chandra. One of the first new studies is the recent one by 
Swartz, Ghosh and Tennant (2003) . Data from this study 
are not yet published, however. 

Finally, we note that another way to probe the nature 
of the objects that constitute ULXs (as well as less luminous 
MBH powered sources) is to look at variability time scales of 



X-ray emission. A tentative linear relation between variabil- 
ity time scale and BH mass has been established although 
with large uncertainties (Markowitz et al. 2003). The rela- 
tion was found to apply to both the SMBHs powering the 
AGN of Seyfert 1 galaxies as well as stellar mass galactic X- 
ray binary systems. Since XRBs display variability across a 
wide range of time scales, this relation, if confirmed, is prob- 
ably only usefully applied to large samples of ULXs. If both 
HMXRBs as well as MBHs power ULXs, variability stud- 
ies might help establish the existence of two distinct time 
scales on which ULXs should be found to display variability 
although with large scatter. 



4.3 Baryonic core accretion at the centre of dwarf 
galaxies and globular clusters ? 

If MBHs do accrete from residual cores of gas or stars, 
then dwarf galaxies and star clusters may be likely places 
to find MBHs and thus luminous X-ray point sources. In 
the hierarchical structure formation scenario, dwarfs orbit- 
ing inside the haloes of larger galaxies correspond to satellite 
subhaloes, and could in principle contain accreting central 
MBHs. There are roughly a dozen known dwarf galaxies 
within the estimated virial radius of the Milky Way, ~ 300 
kpc, for instance. It would be interesting to make a system- 
atic search for unusual X-ray sources at the centres of these 
systems. 
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Closer to the central galaxy, a satellite is likely to have 
been heavily stripped. Previously, we referred to MBHs in 
these systems as 'naked', although they may still be embed- 
ded in a small baryonic core of stars or gas. We would not ex- 
pect to observe anything other than accreting MBHs in these 
systems, unless the core consisted of a significantly number 
of stars bound in a relatively small region. This description 
appears to fit at least some globular clusters (GCs). The ob- 
servations considered above find that ULXs often appear to 
be associated with GCs. Among the formation scenarios that 
are considered for GCs axe that they represent remnants of 
high-redshift pre-galactic fragments (Beasley et al. 2002) 
or even that they were formed in minihaloes that gave rise 
to the first massive stars (Lin & Murray 2002). In the con- 
text of our model, both could correspond to the remnants of 
small satellite systems that formed at high redshift. As such 
they are consistent with the existence of a MBH at their 
centre, accreting from a surviving baryonic core. 

The lack of gas in GCs is not really an obstacle, as we 
already pointed out that a baryonic core of accretable gas 
can in principle be very small. What is more of a problem is 
the required meiss for an MBH to appear as an ULX inside a 
GC. The high-redshift origin of GCs implies that at their for- 
mation MBHs only had a mass equal to the MBH seed mass. 
MBHs of these masses typically only accrete small fractions 
of the their Eddington accretion rate and only in few cases 
this is large enough to power a ULX. However, there is still 
the possibility that since then the MBHs may have grown 
significantly for instance through stellar dynamical processes 
inside the clusters (Ebisuzaki et al. 2001; Mouri & Taniguchi 
2002; Miller & Hamilton 2002; Portegies Zwart & McMillan 
2002). 

4.4 Sources in and around the Milky-Way 

The Milky Way system should correspond roughly to the 
1.6x10^^ Mo halo in our models. For this mass we predict an 
average number of up to three (model C) MBH sources with 
X-ray luminosities exceeding 10^*ergs~^ within the central 
10 per cent of the virial radius, or ~ 30 kpc. A recent survey 
of X-ray sources in the galaxy was carried out by Grimm, 
Gilfanov & Sunyaev (2002) . Their work considers sources 
in the 2.0 - 10 keV band, which gives probably the best rep- 
resentation of the sources; the optical and soft X-ray bands 
are likely to be severely affected by the large hydrogen col- 
umn densities in the disk. The majority of the sources, and 
all of those above lO^^ergs"^, have been positively identi- 
fied with binary systems. The remaining unidentified sources 
have luminosities significantly less than lO'^^ergs"^. Even if 
we assume that the X-ray luminosity of MBHs is close to 
the maximum in the ADAF spectrum, the absence of any 
non-binary ULX in the Milky- Way is still consistent with 
our predictions. On average we would expect about one to 
two ULXs, as can be seen in figure 7, but the standard devi- 
ation is large enough that many haloes would have no such 
sources. 

We can also try to compare results for the entire halo 
of the Milky Way. Figure 10 shows the cumulative X-ray 
luminosity function for the entire halo of a Milky- Way sized 
galaxy. We have chosen a soft X-ray band in the energy range 
0.1 - 2.4 keV to compare with observations of galaxies in the 
Local Group as compiled by Zang & Meurs (2001) . We are 



particularly interested in those Local Group galaxies that 
are satellites of the Milky Way. The observations focus on 
the core regions of these galaxies, which is where we would 
expect MBHs to reside, if they are present at all. For all 
known dwarf satellites of the Milky Way, the X-ray lumi- 
nosity from the cores does not exceed ^ 10'*'^ erg s~^, and 
typically it is much lower. However, for the larger satellites 
the situation is different. For the Large Magellanic Cloud 
(LMC), for example, Kahabka (2002) finds a number of 
sources at X-ray luminosities larger than lO'^^erg s^^ that 
have been classified as binary systems. Since the classifica- 
tion of binaries has been done on the basis of X-ray hard- 
ness ratios this leaves open the possibility that some of them 
could be accreting MBHs. 



5 SUMMARY AND CONCLUSIONS 

We have examined the observational consequences of a pop- 
ulation of free-floating MBHs in galactic haloes. Both our 
predictions and the observations of such sources are highly 
uncertain. In our models, neither the mass nor the abun- 
dance or spatial distribution of MBHs allows for a unique 
identification by a single type of observation alone. As far as 
optical and X-ray signatures are concerned, MBHs could be 
numerous enough to differentiate them from the expected 
count of background sources. On the other hand, even for 
optimistic choices of the model parameters, the mass and 
accretion rates predicted are not large enough for the ac- 
cretion luminosity of MBHs to be of a significantly different 
from that of XRBs. Forthcoming data from Chandra and 
XMM should allow for better spectral modelling of individ- 
ual sources, which may resolve this ambiguity. 

We have shown that our predicted sources, particularly 
from model A, could account for a large fraction of the bright 
emitters inferred from observations. The observed X-ray lu- 
minosity is determined by fitting a spectral model to the 
detected count rate. Depending on the model fitted the re- 
sulting luminosity can differ by a factor of order unity, which 
we have deemed sufficiently accurate to compare with pre- 
dictions. This uncertainty is probably less significant than 
that arising from the assumption made about the interven- 
ing absorbing column densities, which have a strong impact 
on the counts of soft X-ray photons. 

However, while different spectral models may yield sim- 
ilar luminosities they make significantly different assump- 
tions about the underlying emission process and ultimately 
the nature of the ULX source. Typical fitting models include 
a multi-colour disk (MCD) and/or a power law, which have 
the advantage of being specified by few parameters. In most 
cases the colour temperature of the MCD implies a 1 — 10 
M0 BH emitting at or above the Eddington rate. It seems 
that the ADAF model applied to sub-Eddington accreting 
MBHs can match not only the X-ray luminosities but also 
spectral range of detected counts of ULXs. We therefore pro- 
pose to routinely include ADAF models in spectral fits to 
X-ray observations. 

ADAF models do involve one or two more key parame- 
ters, which may introduce degeneracies when used to fit X- 
ray data alone. The advantage is that for a single (as opposed 
to binary) MBH, ADAFs also predict the spectrum in the 
optical regime and as well as the existence of a correlation 




Figure 10. Cumulative X-ray luminosity functions for a 1.6 x 10^^ Mq halo. The shaded area is the standard deviation around the mean 
number given by the lines. Model A is consistent with the absence of MBH sources more luminous than 10^* erg in the Milky- Way 
halo, while for model C this represents a ~ l.5-a deviation. 



of the radio and X-ray fluxes. Consequently, observations in 
these bands could help break the degeneracy arising from 
the increased number of parameters. In contrast for binaries 
these bands may carry the spectral contribution of the visi- 
ble companion, accounting for which adds more parameters, 
too. To fit the overall spectrum ADAFs and composite disk 
models for binary systems thus socm equally well placed. 

Because of the absorption problem optical follow-ups 
arc likely to yield best results for ULXs that are observed 
at high galactic latitudes to avoid the effect of the Galactic 
disk. For such ULXs optical data could help distinguish be- 
tween MBHs and binary systems for sources placed within 
galaxies. To test the prediction of MBHs in the dark mat- 
ter haloes outside galaxies and distinguish them from back- 
ground AGN requires a systematic survey of X-ray point 
sources and optical follow-ups in fields around galaxies. 
Our results in particular are 

• Wc have applied a thin disk and ADAF to model the 
spectral signature of accreting MBHs. Applying this, we pre- 
dict a few ULXs on average per galaxy halo. This is not 
consistent with observations, and therefore indicates either 
that our predicted accretion rates are somewhat too high, 
or that thin disk emission is shut off. 

• B and V band magnitudes reach -10, or -15 when a thin 
disk is included. Because of absorption these are probably 
more useful for follow-ups. 

• For ADAF accretion alone we match the slope of the lu- 
minosity function for one statistical sample of observations, 
and both, slope and normalisation for another one. On aver- 
age we expect to find one ULX in a Milky- Way sized galaxy, 
although, the actual absence of one in the Milky- Way is con- 
sistent with our predictions. 

• Wc did not find any luminosity segregation in the MBH 
sources. In general, MBHs are distributed over a somewhat 
larger range of radii than observed ULXs, but this may be 
partly due to observational selection. 

While many ULXs are clearly linked to star formation 



and young objects such as high-mass XRBs, it is not clear 
that all ULXs need correspond to the same physical sources. 
We have shown that MBHs provide an alternative explana- 
tion for ULXs, and could plausibly account for an appre- 
ciable fraction of observed systems, at least in terms of the 
slope and normalisation of typical X-ray luminosity func- 
tions. MBHs might be particularly interesting candidates to 
explain ULXs in ellipticals, or those seen at large distances 
from star-forming regions in spirals. 

The accretion luminosities predicted here are depend 
sensitively on the mass of MBHs and on the density of ma- 
terial surrounding these systems. As explained in this paper 
and in paper I, we calculate these quantities by assuming 
that seeds merge together efficiently, and that they retain 
around them a core of roughly galactic density. If merging 
does not proceed efficiently and we ignore any significant 
mass increase through gas accretion, then most MBHs will 
have roughly the initial seed mass and X-ray luminosities 
will never reach the ULX range. Similarly, if the debris in 
the core around MBHs had a much lower density then they 
could easily escape detection. 

The luminosities of accreting MBHs in individual bands 
also depend on the parameters of the particular spectral 
model used. For the ADAF model, however, the predicted 
X-ray luminosity is generally proportional to the bolomet- 
ric luminosity. That means we can restate our results more 
generally by saying that if MBHs emit a reasonable fraction 
of their bolometric luminosity in X-rays, they can account 
for many of the observed ULXs. For the ADAF model the 
fraction of the luminosity emitted in X-rays is predicted to 
be fairly high, but in general wc can keep this fraction as a 
free parameter encapsulating the overall uncertainty arising 
from the accretion rate and the particular spectral model 
used. 
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